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Abstract — In this paper, a new filtering method is presents to 
remove Rician noise from magnetic resonance image. This filter 
is based on Neutrosophic set (NS) approach of wiener filtering. A 
Neutrosophic Set (NS), a part of neutrosophy theory, studies 
the origin, nature, and scope of neutralities, as well as their 
interactions with different ideational spectra. Now, we apply 
the neutrosophic set into image domain and define some 
concepts and operators for image denoising. The image is 
transformed into NS domain, described using three membership 
sets: True (T), Indeterminacy (I) and False (F). The entropy of 
the neutrosophic set is defined and employed to measure the 
indeterminacy. The co-wiener filtering operation is used on T and 
F to decrease the set indeterminacy and remove noise. The 
experiments have conducted on simulated MR images from 
Brainweb database with Rician noise added. The visual and the 
diagnostic quality of the denoised image are well preserved. The 
performance of this filter is compared with anisotropic diffusion 
filter (ADF) and unbiased non local mean filter (UNLM). 
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I. Introduction 

Magnetic Resonance Imaging (MRI) is a powerful 

diagnostic technique used in radiology to visualize detailed 
internal structures of the human body. MRI is primarily used 
to demonstrate pathological or other physiological alterations 

of living tissues and is a commonly used form of medical 

imaging [1]. Despite significant improvements in recent years, 
MR images often suffer from low signal to noise ratio (SNR) 

especially in cardiac and brain imaging. The noise in MR 

images follows Rician distribution [2]. This is problematic 
for further tasks such as segmentation of important 
features; classification of images for computer aided 
diagnostics, three dimensional image reconstruction and 
image registration. Therefore, noise reduction techniques are 
of great interest in MR imaging. 


Numerous approaches of denoising MR images have been 
proposed in the literatures [3-12]. In this paper, the 
neutrosophic set approach (NS) of wiener filtering is proposed 
to remove the Rician noise in MRI. First the noisy MRI is 
transformed into the neutrosophic set and the 0) - wiener 
filtering is employed to reduce the indetermination degree of 
the image, which is measured by the entropy of the 
indeterminate subset, after filtering, the noise will be removed. 
The experiments on simulated MR images (Brainweb database 
[13]) having different Rician noise levels demonstrate that this 
approach can perform denoising better. 

II. Neutrosophic MR Image Denoising 

Neutrosophy, a branch of philosophy introduced in [14] as a 
generalization of dialectics, studies the origin, nature and 
scope of neutralities, as well as their interactions with different 
ideational spectra. 

A. Neutrophic Set 

Let U be a Universe of discourse and a neutrosophic set 
A is included in U . An element x in set A is noted as 
x(r, I,F). T,I and F are called the neutrosophic components. 
The element x(7 I. * * * * * 7 ,/,E 7 ) belongs to A in the following way. It 
is t% true in the set, i% indeterminate in the set, and 
/% false in the set, where t varies in T , i varies in I and 
/ varies in F . 

B. Transform the image into neutrosophic set 

Let U be a Universe of discourse and IE is a set of U , 
which is composed by bright pixels. A neutrosophic image 
P NS is characterized by three membership sets T,I,F . a 

pixel P in the image is described as P(T,I,F)md belongs 
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to W in the following way: It is t true in the set, 
i indeterminate in the set, and / false in the set, where 
t varies in T , i varies in 7 and / varies in F . Then the pixel 
P(i 9 j) in the image domain is transformed into the 
neutrosophic set 


domain P NS (; i , j) = {T (; i , j), I(i, j), F(i, j)} . T (i, j), I(i, j ) 
and F(i, j) are the probabilities belong to white pixels set, 
indeterminate set and non white pixels set respectively, which 
are defined as: 
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where g(i 9 j) is the local mean value of the pixels of the 
window. S(i,j) is the absolute value of difference between 
intensity g(i 9 j) and its local mean value g(i 9 j) . 

C. Neutrosophic image entropy 

For a gray image, the entropy is utilized to evaluate the 
distribution of the gray levels. If the entropy is the maximum, 
the intensities have equal probability. If the entropy is small, 
the intensity distribution is non-uniform. 

Neutrosophic entropy of an image is defined as the 
summation of the entropies of three subsets T , 7 and F : 
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where En T , Enj and En F are 

the entropies of sets 


7 1 , 7 and F respectively. p T (/) , pj (/) and p F (/) are the 
probabilities of elements in T 7 , 7 and F respectively, whose 
values equal to i . 


D. CO -wiener filtering operation 

The values of 7(7, j) is employed to measure the 
indeterminate degree of element 7^ (i 9 j) . To make the set 
7 correlated with T and F , the changes in T and F 
influence the distribution of element in 7 and vary the entropy 
of 7 . A co - wiener filtering operation for P NS , P NS (a?) , is 
defined as: 
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where Sj, (/, j) is the absolute value of difference between 

intensity f (/, j) and its local mean value T (/, /) at 
(/, /) after & - wiener filtering operation. 

The MRI denoising method based on neutrosophic set 
approach of wiener filtering is summarized as below: 

Step 1 : Transform the image into NS domain; 

Step 2: Use co - wiener filtering operation on the true subset 
T to obtain T w ; 

Step 3 : Compute the entropy of the indeterminate subset 

ico,En } ( 0 ; 

1 co 

En ~ [ ( i + X) — En- (i) 

Step 4: if — — < 8 , go to Step 5; 

En- O') 

1 (0 

Else T — f M , go to Step 2; 

Step 5: Transform subset T (0 from the neutrosophic domain 
into the gray level domain. 



III. Results and Discussion 


The experiments are conducted on simulated MR images 
obtained from the Brainweb database [13]. The dataset 
consists of T1 and T2- weighted MR images. The experiments 
are performed on various images (T1 and T2- weighted 
normal MR images; T1 and T2- weighted MR images with 
Multiple Sclerosis) degraded by different noise levels. Here, 
the performance of the proposed method is demonstrated for 
T1 -weighted MR image of normal brain with 9% noise level 
and T1 -weighted brain MR image with MS lesion corrupted 
by 9% noise level. This method is compared with Anisotropic 
Diffusion Filter (ADF) [3] and Unbiased Non Local Mean 
(UNLM) filter [11]. The results of the different denoising 
filters on a T1 -weighted MR image of normal brain and with 
MS lesion corrupted by 9 % Rician noise level are shown in 
Fig. 1 and Fig. 2 respectively. From these results, the proposed 
denoising filter produce good denoised image in terms of 
visual perception. 



(e) 


Figure 1. Results of different denoising filters on a T1 -weighted MR image 
of normal brain with 9% noise level. From top to bottom: (a) the original 
image, (b) the noise corrupted image, denoised images using (c) AD filter, (d) 
UNLM filter and (e) NS wiener filter. 



Figure 2. Results of different denoising filters on a T1 -weighted MR image 
of MS lesion brain with 9% noise level. From top to bottom: (a) the original 
image, (b) the noise corrupted image, denoised images using (c) AD filter, (d) 
UNLM filter and (e) NS wiener filter 


The performance of the denoising algorithm is measured by 
using the quality metrics such as peak-signal-to-noise ratio 
(PSNR), the Structural Similarity (SSIM) index [15] and the 
Quality Index based on Local Variance (QILV) [16]. The peak 
signal to noise ratio in decibel (dB) is measured using the 
following formula: 


PSNR = -10 


log 


Z i—H- 
i = 0 


H xW x255 2 


(18) 


where j) and I d (i, j ) represent the intensities of pixels 
(/. /) in the original image and denoised image respectively. 


The higher the PSNR, the better the denoising algorithm is. 
SSIM and QILV give the measure of the structural similarity 
between the original and the denoised images and are in the 
range of 0 to 1. The SSIM works as follows: Let xand y be 













two non negative images, where as one has perfect quality. 
Then, the SSIM can serve as a quantitative measure of the 
similarity of the second image. The system separates the task 
of similarity measurement into three comparisons: luminance, 
contrast and structure. It can be defined as 
[l/J/J + Q )(2cr + C 2 ) 
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a are the standard deviations respectively, o can be 


estimated as 


xy 




( 20 ) 


C x and C 2 in eqn. 19 are constants and the values are given as 

C x =(K X L) 2 and C 2 ={K 2 L) 2 where K X ,K 2 « 1 is a small 
constant and L is the dynamic range of the pixel values (255 
for 8 bit gray scale images). 

The QILV between two images x and y can be defined as 


QILV(x,y) = 
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where ju Vx and fi Vy are the estimated means of the local 
variance and G Vx and G Vy are the standard deviations of the 


local variance respectively. G VxVy 
the two images. 


is the covariance between 


The efficiency of the proposed denoising filter is compared 
with ADF and UNLM. The experiments carried on a simulated 
T1 -weighted MR image of normal brain with 1%, 3%, 5%, 
7%, 9%, 15% and 21% of Rician noise levels are added. The 
quality metrics PSNR, SSIM and QILV are calculated and the 
comparison of these metrics for ADF, UNLM and proposed 
filter is given in Fig. 3 and Fig. 4. The PSNR, SSIM and QILV 
values obtained for the T1 and T2 weighted images using the 
aforementioned denoising techniques are tabulated in Table I. 
From these results, the proposed filter provides better denoising 
by obtaining high PSNR, High SSIM and high QILV for all the 
values of noise levels. 



(a) 



(b) 



TABLE I 

Comparison of the Denoising Techniques Based on the Performance 
Metrics for Different MR Images 


MR Image 

Denoising 

Methods 

Performance Metrics 

PSNR(dB) 

SSIM 

QILV 

Tl - Weighted brain MRI 
corrupted by 7% Rician 
noise 

ADF 

23.98 

0.962 

0.995 

UNLM 

29.72 

0.9908 

0.7309 

NS Wiener 

31.91 

0.9943 

0.9871 

Tl - Weighted brain MRI 
with MS lesion corrupted 
by 9% Rician noise 

ADF 

21.59 

0.9331 

0.7702 

UNLM 

24.59 

0.9714 

0.9692 

NS Wiener 

29.95 

0.9903 

0.993 

T2- Weighted brain MRI 
corrupted by 15% Rician 
noise 

ADF 

18.24 

0.8808 

0.6555 

UNLM 

23.44 

0.9532 

0.9441 

NS Wiener 

24.08 

0.9682 

0.9882 


(c) 

Figure 3. Comparison of the denoising results obtained from ADF, UNLM, 
and NS Wiener Filters based on the quality measures for Tl- weighted normal 
brain MRI: (a) PSNR (b) SSIM (c) QILV. 

IV. Conclusion 

In this article, a novel MRI denoising technique based on 
neutrosophic set approach of wiener filtering has been 
proposed. The image is described as a NS set using three 
membership sets T, I and F. The entropy in neutrosophic image 
domain is defined and employed to measure the 
indetermination. The wiener filter is applied to reduce the set’s 
indetermination and remove the noise in the MR image. The 
performance of the proposed denoising filter is compared with 
ADF and UNLM based on PSNR, SSIM and QILV. 
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Figure 4. Comparison of the denoising results obtained from ADF, UNLM, 
and NS Wiener Filters based on the quality measures for Tl- weighted brain 
MRI with MS lesion: (a) PSNR (b) SSIM (c) QILV 


The experimental results demonstrate that the proposed 
approach can remove noise automatically and effectively. 
This filtering method tends to produce good denoised image 
not only in terms of visual perception but also in terms of the 
quality metrics. 
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